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I.  Purpose 


Computer  codes  modeling  physical  phenomena  may  be  divided  into 


three  categories  as  depicted  in  Figure  1-1. 


Fig.  1-1  Computer  Code  Models 


The  first  category  is  the  full  physics  model  which  contains  the 
most  rigorous  scientific  development,  adhering  to  generally  accepted 
theory  in  every  possible  detail.  The  second  category  is  a  model  based 
on  physical  principles  but  containing  approximations.  This  type  of 
model  is  generally  less  sophisticated  than  the  full  physics  model  and  is 
easier  and  les3  costly  to  run.  The  third  category  is  the  model  based  on 
empirical  fits  of  data.  An  empirical  model  is  very  fast  running  and 
inexpensive  to  use. 

The  purpose  of  this  project  is  to  produce  a  set  of  empirical 
equations  to  model  the  fractional  rate  of  arrival  of  radioactivity  on 


the  ground.  The  model  will  be  limited  to  fallout  from  n  nuclear  burst 
at  the  surface  of  the  eartfi  or  close  enough  such  that  the  fireball 
touches  the  ground. 


II .  Background 


Since  the  commencement  of  peacetime  atmospheric  testing  of  nuclear 
weapons,  residual  radiation  has  been  recognized  as  a  hazardous  nuclear 
effect.  Even  before  the  first  detonations  which  produced  significant 
local  fallout,  the  concept  of  fallout  was  known  and  published  (11:15). 
In  the  intervening  years  numerous  models  have  been  developed  to  predict 
radiation  fallout  characteristics.  The  models  may  be  classified  as 
either  disk  tossers  or  smear  models. 

Disk  Tossers.  The  disk  tosser  fallout  prediction  model  discretizes  the 
nuclear  debris  cloud  into  elements  (at  least  one  dimension  of  particle 
size  and  one  of  space  (3:205))  and  follows  the  time  history  of  each 
element.  Superposition  of  all  grounded  elements  yields  the  fallout 
footprint.  Because  of  the  numerical  quadrature  of  disc  tosser  codes, 
computers  are  required  to  complete  the  calculations.  The  Department  of 
Defense  Land  Fallout  Prediction  System  (Ref.  5),  DELFIC,  is  a  disc 
tosser  code  and  is  generally  regarded  as  the  foremost  of  such  codes.  It 
is  a  full  physics  model  and  uses  a  quadrature  of  up  to  10,000  elements 
to  describe  the  local  fallout  pattern  (7:27).  Such  required  fine  meshes 
are  expensive  in  computer  time  and  money. 

Smearing  Codes.  Smearing  codes  do  not  discretize  the  stabilized  debris 
cloud  into  a  fine  mesh  for  numerical  transport  of  each  element  until 
ground  impact,  but  rather  continuously  deposit  the  activity  in  the 
stabilized  cloud  along  the  ground  until  all  radioactivity  of  interest  is 
grounded.  The  process  may  be  compared  to  smearing  a  piece  of  chalk 
continuous ly  a  long  a  chalkboard  until  the  initial  piece  is  entirely 


depleted.  Hence  the  label  of  smearing  codec.  "igure  IT-1  is  an 
illustration  of  the  smearing  process. 


Figure  XI —1  Smearing  a  Radioactive  Debris  Cloud 


One  of  the  earliest  smearing  codes  was  created  by  the  Weapons 
Systems  Evaluation  Group  in  the  Pentagon.  The  heart  of  tin;  code,  called 
WSEG-10,  lies  in  a  function  g(t)  (10:208),  which  is  defined  in  the  WSEG- 
10  report  as  follows  "Specifically,  for  any  radioactive  cloud 
configuration  there  exists  some  function  of  time,  g(t),  that  represents 
the  fraction  of  total  radioactivity  that  arrives  on  the  ground  per  unit 
time."  (9:4).  The  smearing  equation  using  g(t)  is 

t 

A(x,y,t)  =  AT ( t )/  f( x ,y , t ' ) g( t ' ) d t *  [Ref  d  eqn  (?)]  (l) 

o 

where  A(x,y,t)  is  the  activity  footprint  on  the  ground 
in  Curios/m' 


I  I 


A,p(t)  is  the  total  radioactivity  in  the  cloud  at  time  t 
in  Curies 

f(x,y,t')  is  the  normalized  distribution  of  activity  per 
area  horizontally  across  the  cloud  at  time  t 
in  1/m^ 

The  integral  term  smears  the  horizontally  distributed  activity  cloud 
along  the  ground  at  a  rate  determined  by  g(t)  (3:208). 

WSEG-10  has  been  in  use  continuously  for  the  last  20  years.  Its 
popularity  is  due  partly  to  its  simplicity  and  fast  computation  time 
(2:111-2).  WSEG-10  uses  a3  the  g( t)  function  a  simple  decreasing 
exponential  which  was  fit  to  nuclear  test  fallout  data. 

But  WSEG-10  has  come  under  strong  criticism  by  Russell  (10:209)  for 
its  treatment  of  particle  size  distribution,  fractionation,  and  g(t);  by 
Polan  (8:31)  for  its  [lack  of]  treatment  of  particle  .sizes  and  settling 
rates;  and  by  Tlorment  (7:103)  for  the  slope  of  the  particle  size  - 
activity  distribution.  Fractionation  is  the  process  whereby  part  of  the 
radioactivity  in  the  stabilized  debris  cloud  is  1  i.-st  riba  ted  throughout 
the  volume  of  fallout  particles  .ami  the  remainder  of  the  radioactivity 
is  distributed  on  the  surface  of  the  particles.  Bridgman  and  Bigelow 
have  published  a  new  fallout  prediction  smearing  model  'Ref.  0  which 
derives  the  g(t)  function  from  physics  and  a  result  is  dependent  on 
particle  size  -  activity,  particle  settling  rut«s  an i  accounts  for 
fractionation  (3:2r  .) 


The  Bridgman  -  1 


jdel  is  the  model  is<*d  n  this  project. 


III.  Method 


$ 


Theo  ry.  If  a  particle  ac ti v i ty- a i ze  distribution,  A(r),  can  be 

determined,  the  fraction  of  activity  arriving  on  the  ground  must  be 

related  to  A(r)  by  (3:21 0) 

g(  t)  dt  =  -A( r)  dr  (?) 

The  minus  sign  indicates  particles  of  decreasing  radius  are  reaching  the 

ground  as  time  increases.  Then 

g(t)  =  -A(r)  dr/dt  (3) 

For  a  particle  number-size  distribution,  N(r),  which  is  taken  to  be  a 

log-normal  distribution,  calculation  of  the  activity-size  distribution 

is  straight  forward.  The  log-normal  distribution  is 

N(r)  =  NT//(2n)Br  exp(-1 /?)[ (l n( r)-aQ )/p]  (4) 

where  aQ  is  the  logarithm  of  the  median  particle  radius 
3  is  the  geometric  deviation 
Nrp  is  the  total  number  of  particles 

If  all  of  the  activity  in  the  stabilized  cloud  were  distributed  in  the 
volume  of  the  fallout  particles  then  the  activity-size  distribution 
would  be  proportional  to  the  third  moment  of  N(r).  If  all  of  the 
activity  were  distributed  only  on  the  surface  of  the  particles  then  A(r) 
would  be  proportional  to  the  second  moment  of  N(r).  A  property  of  log¬ 
normal  distributions  i3  that  the  nth  moment  of  such  a  distribution  is 
itself  a  log-normal  distribution  with  the  same  geometric  deviation,  B, 
and  an  a  defined  by  the  relation  (1:12) 

an  =  a0+  nB  (3) 

where  ac  is  the  logarithm  of  the  median  particle  radius.  Since  a 
portion  of  the  radioactivity  in  the  stabilized  cloud  is  volume 
distributed  and  the  remainder  is  surface  distributed ,  the  activity-3ize 
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distribution  can  be  expressed  as  a  weighted  sum  of  two  log-normal 
distributions 

A( r)  =  FvAT/[/l2n)Br]  exp(-1 /2)[ (l n( r)-Oj)/3]2 

+  (1  -Fv)AT/[/(2ir)  Br]  exp(-1 /2) [  ( 1  n(  r)-a2 )/6]2  (6) 

where  A<p  is  the  total  activity  at  time  t 

Fv  is  the  fraction  of  activity  that  is  volume  distributed 

Using  appropriate  fall  mechanics  for  particles  descending  in 
atmosphere  the  time  for  a  particle  of  radius  r  to  reach  the  ground  from 
an  initial  altitude  can  be  determined,  and  dr/dt  can  be  found. 

Model.  A  fast  running  computer  code  based  on  the  smearing  model 
developed  by  Bridgman  and  Bigelow  (Ref.  3)  was  used  to  generate  g(t) 
data  for  this  project.  The  Bridgman-Bigelow  model,  hereafter  called  the 
AFIT  smear  model,  was  chosen  because  it  accounts  for  fractionation  and 
it  contains  a  g(t)  function  which  was  derived  from  physical  principles 
rather  than  empirical  fits  to  available  test  shot  data.  Further,  the 
g(t)  values  calculated  with  this  model  have  been  shown  to  agree  well 
with  DELFIC  predictions  (1 sVI  — 2).  Figures  1 1 1  —  1  through  III-5 
illustrate  this  agreement. 

The  Bridgman  -  Bigelow  method  was  recoded  by  this  author  and  the 
source  program  is  listed  in  Appendix  A.  Inputs  to  this  code  are 
particle  size  distribution  parameters,  f  rac  t  iona  t.ion  ratio  and  the 
desired  yield.  The  program  will  compute  g(t)  data  for  fallout  to  any 
time  after  cloud  stabilization.  ^or  purposes  of  this  project  data  was 
limited  to  local  fallout  where  local  fallout  is  defined  as  that 
radioactivity  which  reaches  the  ground  within  24  hours. 
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Fig.  III-4  Comparison  of  q(t) 
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Procedure.  For  each  particle  size  distribution  listein  Table 
III— 1 »  g(t)  data  were  generated  for  five  different  nuclear  yields.  All 
of  the  distributions  are  log-normal.  The  five  yields  chosen  were  1 
Kiloton  (KT),  10KT,  100KT,  1  Megaton  (MT)  and  15MT.  The  four  smaller 
yields  are  standard  yields  of  interest.  The  15MT  yield  was  chosen 
because  of  the  fallout  data  available  from  test  shot  Castle  Bravo  where 
the  yield  was  approximately  15MT  ( 4:37 )•  The  generated  g(  t)  data  were 
then  fit  with  a  Laurent  series  for  each  yield 


g(t)  =  l  Ci  t1-i  (7) 

i=1 

where  C.^  is  the  coefficient  of  the  ith  term 

t  is  the  time  of  activity  arrival  on  the  ground 


Table  TII-1 

Particle  Size  Distributions 


Distribution 

Name 

Median 

Radius 

(pm) 

Geometric 

Deviation 

Fractionation 

Value 

DELFIC  Default 

.204 

4.0 

.68 

High  Yield 

Standard 

1 1  .4 

2.16 

.74 

Low  Yield  Standard 

.157 

4.0 

.  66 

NRDL-C61 

.0424 

5.576 

.74 

Geometric  Deviation  and  Fractionation  Value  are  dimensionless. 

The  DELFIC  default  distribution  is  the  particle  size  distribution  used 
by  the  DELFIC  disc  tosser  code  on  the  absence  of  other  operator  input. 
The  remaining  distributions  were  used  by  Norment  in  a  PHI. FTC  sensitivity 
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study  (6:25)*  The  High  Yield  Standard  distribution  was  intended  to 
represent  the  distribution  used  by  W  S EG  —  10*  The  Low  Yield  S t n d r d 
distribution  was  chosen  to  represent  the  particle  size-activity  data 
from  the  Small  Boy  test  shot  (2:11-13)*  The  NRDL-Cbl  distribution  was 
used  by  the  Navy  Radiological  Defense  Laboratory  to  describe  South 
Pacific  Coral.  The  fractionation  values  are  taken  from  work  by  Bigelow 
(2:IV-3). 

Table  III-2  is  a  compilation  of  the  Ci's  for  the  Laurent  fit  to  the 
AFIT  smear  code  g(t)  data  generated  by  using  the  DELPIC  default  particle 
size  distribution.  A  polynomial  least  squares  fit  was  made  of  each  set 
of  Laurent  coefficients,  that  is,  the  five  Cl's,  the  five  C2' s,  etc.  as 
follows : 

Ci  -  l  Kj  [ ln(YM) ] (3) 

j  =  1 

where  C^  is  the  Laurent  coefficient  being  fit 
K-  i3  the  polynomial  coefficient 
Ym  is  the  nuclear  yield  in  megatons 
n  is  the  degree  of  polynomial  required  to 
fit  the  Laurent  coefficient 

Table  III-3  is  a  matrix  of  K-'s  for  the  DELFIC  default  distribution 
Laurent  coefficients.  Appendix  B  contains  tabular  listings  of  the 
Laurent  and  polynomial  coefficients  for  all  the  size  distributions 
analyzed. 

The  benefit  of  fitting  the  Laurent  coefficients  with  a  polynomial 
is  that  g(t)  data  can  then  be  calculated  for  any  arbitrary  nuclear  yield 
between  1KT  and  1  5  MT  for  any  of  the  four  particle  size  distributions. 
An  illustrative  example  follows. 

Find  the  fractional  activity  deposition  rate,  g(t),  at  the  arrival 


time  of  2  hours  given  a  50KT  yield  (.OSM?)  and  using  the  OKI, FTC  default 


particle  size  distribution.  Using  equation  3, 

Cl  =  K1*[ln(.05)]°  +  K2*[ln(.05)]1  +  K3*[ 1 n( .05 ) ]2 
+  K4*[ln(.05)]3  +  K5[ln(.05)]4 
Cl  =  -.001 17313*(-2.996)°  +.00004020*(-2. 996)1 
+  .00000601  *(-2.996 )2  -.00000427*( -2.996)3 
-.00000077* (-2. 996) 4 

Cl  =  -.00117313  -.00012044  +.00005395  +.00011483  -.00006204 
Cl  =  -.00118683 

C2  through  C7  are  calculated  in  the  same  manner.  Using  equation 
the  calculated  Laurent  coefficients 

g(2  hr)  =  -.00118683/1  +.14255015/2  +.08281579/4  -.08754622/8 
+.03127940/16  -.00532854/32  +.00035457/64  per  hour 
g( 2  hr)  =  .0817  per  hour 
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IV.  Results 


Accuracy  of  Laurent  Series  Fit.  Four  particle  si.se  distributions  were 
used  to  generate  g(t)  data  for  five  different  nuclear  yields.  The  g( t) 
data  for  each  yield  in  each  distribution  were  fit  with  a  Laurent  series. 
To  achieve  a  fit  of  good  accuracy,  g(  t)  data  from  .2  hours  to  24  hours 
after  cloud  stabilization  time  were  used,  and  coefficients  for  a  sixth 
degree  Laurent  series  were  computed.  3y  ignoring  data  from  times 
earlier  than  .2  hours  the  Laurent  series  coefficients  will  produce  g(t) 
data  that  is  within  approximately  4  percent  of  the  smearing  code  g(  t) 
data  in  most  cases.  Ignoring  the  data  from  times  earlier  than  .2  hours 
is  an  insignificant  restriction  because  radioactive  fallout  is  rarely  a 
concern  so  soon  after  a  nuclear  detonation.  Other  nuclear  effects  are 
of  greater  and  overriding  concern. 

The  only  exception  to  the  4 %  agreement  over  the  range  .?  <  t  <  24 
hours  is  the  1 KT  yield  of  the  High  Yield  Standard  particle  size 
distribution.  The  High  Yield  Standard  distribution  produces  a  g(  t) 
function  characterized  by  a  very  sharp  peak  followed  by  a  slowly 
dropping  plateau.  The  achieve  a  Laurent  fit  with  accuracy  of  order  4 
percent  it  was  necessary  to  fit  two  different  time  intervals.  The  time 
intervals  chosen  were  from  .3  to  2  hours  and  f rom  ,?  to  24  hours.  Table 
I V  —  1  illustrates  the  accuracy  of  the  Laurent  fit  for  each  of  the  five 
yields  in  each  distribution.  The  g( t)  data  generated  using  the  Laurent 
fit  are  compared  to  the  g(t)  data  computed  with  the  AR1T  smearing  model. 
Tabular  lists  of  the  Laurent  scries  coefficients  are  provided  in 
Appendix  3. 


Table  IV-1 

Laurent  Fit  to  Smear  g(t)  Data 


Distribution  Name 

Yield 

Percent  Error* 

DELFIC  Default 

15MT 

2.62 

1 MT 

2.65 

100KT 

5.02 

10KT 

2.07 

1KT 

3.14 

High  Yield  Standard 
•3-2  Hours 

15MT 

1  .06 

1  MT 

1.27 

100KT 

0.88 

10KT 

1  .03 

1KT 

0.20 

2-24  Hours 

15MT 

1  .50 

1  MT 

0.70 

100KT 

2.56 

1 0KT 

4.02 

1  KT 

53-0 

Low  Yield  Standard 

15MT 

1 .05 

1  MT 

0.83 

100KT 

2.67 

1 0KT 

3.60 

1  KT 

2.28 

NRDL-C61 

15MT 

1.25 

1MT 

0.55 

100KT 

1.23 

10KT 

2.92 

1  KT 

2.28 

*  Absolute 

error  is  noted.  Actual  percent 

error  is 

sometimes 

positive  or  negative  due  to  the 

oscillatory 

behavior 

of  a  high  degree  Laurent  series 

fit. 

Accuracy  o f  Polynomial  Fit  to  Lauren t  Coefficients.  Each  set  of  Laurent 
series  coefficients  for  the  five  chosen  yields  in  each  particle  size 
distribution  was  fit  with  a  polynomial  series  as  a  function  of  yield. 
That  is,  the  set  of  five  Cl's,  five  C2's,  etc  for  each  distribution  were 
fit  with  polynomials.  The  polynomial  fits  to  the  Laurent  coefficients 
provide  the  capability  to  generate  g(t)  data  for  any  nuclear  yield 
between  1KT  and  15MT.  To  achieve  an  accurate  fit  coefficients  for 
polynomials  to  4th  degree  were  computed.  Coefficients  for  3rd  degree 
polynomials  were  computed  for  all  four  particle  size  distributions  and 
coefficients  for  2nd  degree  polynomials  were  computed  for  the  DELFIC 
default  and  NRDL-C61  distributions.  In  all  cases  the  higher  order 
polynomial  yields  the  most  accurate  computation  of  g^t)  data,  although 
satisfactory  results  may  be  obtained  from  second  order  polynomials  in 
some  cases.  Table  IV-2  is  a  listing  of  the  accuracy  of  g(t)  data 
calculated  with  the  polynomial  fits  to  the  Laurent  coefficients  as 
compared  to  the  g(t)  data  generated  by  the  AFIT  smear  model.  Tabular 
lists  of  the  polynomial  coefficients  are  provided  in  Appendix  B. 
Coefficients  for  2nd  degree  polynomials  are  listed  for  those 
distributions  where  the  accuracy  of  the  g(  t)  data  for  the  3rd  degree 
polynomial  is  generally  within  percent  of  the  AFIT  smear  model  data. 

In  Appendix  B  Figures  B-1  through  B-20  graphically  compare  the  g(t) 
data  for  each  of  the  five  yields  for  each  particle  size  distribution. 
Each  figure  displays  the  g( t)  data  a3  computed  by  the  AFIT  smear  model, 
the  Laurent  fit  to  the  smear  model  and  all  applicable  polynomial  fits  to 


the  Laurent  coefficients. 


Table  1V-2 

Accuracy  of  Polynomial  Pit  g( t)  Data 


Distribution  Name 


DELFIC  Default 


High  Yield  Standard 
•  3-2  Hours 


2-24  Hours 


Low  Yield  Standard 


%  Error  %  Error 
4th  Deg.  3rd  Deg. 
Polynom.  Polynom. 


Yield 


15MT 
1 MT 
100KT 
1 0KT 
1KT 


15MT 
1  MT 
100KT 
1 0KT 
1KT 


15MT 
1  MT 
100KT 
1 OKT 
1  KT 


15MT 
1HT 
100KT 
1 0KT 
1  KT 


15MT 
1  MT 
1 0OK'T 
1 OKT 


2.62 

2.64 

3-02 

2.07 

3.09 


1.13 
1 .00 
0.33 
0.25 

0.21 


1.49 

0.70 

2.56 

4.06 

54.41 


2.57 

0.83 

2.67 

3.60 

3.24 


1 .24 

0.55 

1.22 

2.87 


2.45 
3-48 
2.82 
1 .51 

2.80 


5.71 

12.13 

9.35 

2.97 

0.73 


1 .67 
1  .37 
1  .67 

6.56 

49.60 


7.11 
2.14 
4.  SI 

33-21 

135-7 


1 .30 

0.94 

1  .52 

2.71 


%  Error 
2nd  Deg 
Polynom 


3.49 

4.49 
2.62 
5.00 
3.58 


N/A 


N/A 


N/A 


2.16 

4-30 

1.21 

5.67 


NRDL-C6 1 


V.  Conclusions  and  Recommends t i ons 

Conclusions.  The  fractional  rate  of  radioactivity  deposition,  g(t),  can 
be  calculated  using  an  empirically  derived  Laurent  series  equation.  A 
polynomial  function  of  nuclear  yield  can  be  used  to  calculate  Laurent 
coefficients  which  are  in  turn  used  to  calculate  g(t).  Values  of  g(t) 
thus  calculated  will  be  within  4  percent  of  values  produced  by  the  APIT 
smear  model.  An  exception  is  calculation  of  low  yield  g(  t)  data  using 
the  High  Yield  Standard  particle  size  distribution.  Calcilation  of  g(t) 
data  can  be  made  at  an  accuracy  of  4  percent  for  nuclear  yields  between 
1  KT  and  15MT.  This  range  i3  sufficiently  broad  for  aceommodat ion  of 
tactical  and  strategic  applications. 

Use  of  4th  degree  polynomials  will  produce  g(  t)  values  which  most 
closely  approach  the  values  produced  by  the  APIT  smear  model.  Use  of 
lower  degree  polynomials  yields  percentage  errors  of  less  than  5  percent 
in  some  cases,  as  shown  in  Table  IV-?. 

No  attempt  wa3  made  to  extend  the  useful  range  beyond  the  limits  of 
1  KT  or  15MT.  In  general,  attempts  to  calculate  g(  t)  data  for  yields 
outside  this  range  may  result  in  large  percentage  errors. 

Values  of  g( t)  may  be  obtained  with  the  aid  of  any  hind  calculator. 
Large  computer  systems  like  those  required  for  DPI. PM  or  small  systems 
like  those  required  for  fast  running  smearing  codes  such  as  the  APIT 
smear  model  are  not  required.  The  capability  to  determine  fractional 
fallout  rates  without  computer  assistance  is  advantageous  to  military 
field  commanders  and  planners. 

Recommendations.  Only  four  particle  size  distributions  were  analyzed. 
These  distributions  were  developed  to  represent  Nevada  test  site  soil 
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and  South  Pacific  coral.  Other  particle  size  distributions  for  soils 
which  are  of  interest  to  military  commanders  and  planners  should  be 
developed  and  analyzed. 
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This  appendix  contains  the  computer  listing  of  the  code  "PROGRAM 
g(t)".  It  was  written  in  BASIC  for  a  Nippon  Electric  Corporation  (NEC) 
personal  computer. 

The  code  was  written  to  generate  the  AEIT  smear  g(t)  data  for  this 
project.  Requirements  were  to  provide  g(t)  data  for  five  specific 
yields.  Using  the  pancake  approximation  of  the  nuclear  debris  cloud 
(3:214)  would  necessitate  reading  into  the  code  not  more  than  10  sets  of 
Laurent  coefficients.  In  the  interest  of  flexibility  and  generality, 
all  the  Laurent  coefficients  were  read  into  the  program  which  are 
required  to  enable  calculating  the  g(  t)  data  for  any  yields  up  to  15MT. 

The  following  is  a  glossary  of  variables  used  in  the  code. 


Variables 


AO  -  The  log  of  the  mean  particle  radius 

A2  -  Log-normal  distribution  parameter  of  the  second  moment 
of  the  particle  size  distribution 

A3  -  Log-normal  distribution  parameter  of  the  th  i.  rd  moment 
of  the  particle  size  distribution 

AR  -  Activity  of  particle  of  radius  R 

B  -  Log  of  the  geometric  deviation  of  the  particle  size 
distribution 

Cl  through  C7  -  Laurent  coefficients  for  R  =  r(t,Z)  for 
AFIT  smear  model  (3:214) 

DR  -  rate  of  change  of  radius  of  particle  reaching  the 
ground  in  interval  dt 

FV  -  Fractionation  ratio 


PROGRAM  g(t) 

This  program  is  a  smearing  fallout  prediction  code  for 
the  calculation  of  g(t)  data.  The  code  is  baaed  on  the 
APIT  smearing  model.  The  program,  as  written,  will 
generate  g(t)  values  for  five  different  nuclear  yields. 
The  yields  are  1 KT,  10KT,  100KT,  1 MT,  and  15MT. 

The  particle  size  distribution  parameters  must  be  input. 

INPUT  "DISTRIBUTION  NAME"  ;NS 

INPUT  "MEAN  PARTICLE  RADIUS,  [um]"  ;  PR 

AO  =  LOG(PR) 

INPUT  "GEOMETRIC  DEVIATION,  0"  ;GD 
B  -  LOG(GD) 

INPUT  "FRACTIONATION  RATIO  "  ;  FV 


PI  =  3-14159 


DIM  Cl 
DATA  - 

DATA  - 

DATA  - 

DATA  - 

DATA  - 

DATA  - 


DATA  - 


DATA  - 


DATA  - 


DATA  - 


DATA  - 


DATA  - 


(85)  ,C2(85) 

•42494E-1 5, 

.831 23E-06 , - 

•13323E-13, 

.16550E-05,- 

•99099E-13, 

•24713E-05,- 

•40897E-12, 

•32801 E-05,  - 

• 1 2220E-1 1 , 

.4081 2E-05 , - 

.29772E-1 1 , 

•48748E-05,- 

•  63023E-1 1  , 

•56610E-05,- 

.12036E-10, 

•64397E-05,- 

.21 240E-1 0, 

.7211  IE-05, - 

•35220E-10, 

•79750E-05,- 

•55534E-10, 

.8731 7E-05 , - 

.84005E-10, 

•94814E-05,- 

.1227 IE-09, 

.  1Q224E-04,- 
. 1 7400E-09, 
.10959E-04,- 
.24045E-09, 

. 1 1 688E-04 , - 
•32495E-09, 

. 1 2 409 E -04 , - 
.43064E-09, 


,C3(85)  ,C4 

37874E-1 2, - 

17840E-05, 

59598E-1 1 , - 

17823E-05, 

29668E-10,- 

17806E-05, 

92185E-10,- 

17789E-05, 

22123S-09,- 

1777 IE-05, 

45092E-09,- 

17754E-05, 

82136E-09,- 

17744E-05, 

1 3778E-08, - 
17740E-05, 
21 698E-08, - 
17732E-05, 
32509E-08,- 
17721 E-05, 
46784E-08,- 
17709E-05, 
65129E-08,- 
17695E-05, 
881 7  IE-08, - 
17681 E-05, 

1 1 656E-07 , - 
17664E-05, 

1 5094E-07, - 
17647E-05, 

1 9202E-07, - 

17630E-05, 

24047E-07,- 


(85)  ,C5(85) 

.99764E-10, 

.131 10E-04 

.78875E-09, 

.18523E-04 

.26305E-08, 

.22664E-04 

.61 605E-08, 

.26144E-04 

.1 1887E-07, 

.29202E-04 

. 20291 E-07, 

.31 958E-04 

.31838E-07, 

.34484E-04 

•46963E-07, 

.36829E-04 

.66068E-07, 

•39024E-04 

.89534E-07, 

.41 093E-04 

. 1 1772E-06, 

.43055B-04 

. 1 5098E-06 , 

.44924E-04 

.18962E-06, 

.4671  IE-04 

•23594E-06, 

.48425E-04 

.28419E-06, 

•50073E-04 

•34063E-06, 

. 51 662E-04 
•40353E-06. 


,C6(S5)  ,07(85) 
•  23375E-07 , 

•92734E-07, 

.20692E-06, 

.36479E-06, 

. 5651 6E-06 , 

.80695E-06, 

.  1 0892E-05 , 

.  14107E-05, 

.  1 7705 E-05 , 

.21673E-05, 

. 26002E-05 , 

. 30682E-O5 . 

.36702E-05, 

.41 052E-05 , 

.46722E-05, 

.52701 E-05, 

•58983E-05, 


. 1 31 23E-04,-. 1761  IE-05,  .531 963-04 
370  DATA  -.560873-09,  . 29700B-07 , - •  473085-06 ,  . 6y,y>:-91i 
. 1 3831 E-04,-. 1 7592E-05 ,  .54681 S-04 
380  DATA  -.719235-09,  . 36229B-07 ,  - . 549475-06 ,  .774  1  1  5-05 
.14531 E-04 , - . 1 7572E-05 ,  . 56 1 21 E-04 
390  DATA  -.90951E-09,  . 43702E-07 ,  - . 63206E-06 ,  .7l>‘-3  .'E-05 
. 1 5225E-04 ,  - .  1 7552E-05 ,  • 575 1 8E-04 
400  DATA  -.11 358E-08,  .521 895-07, -.72345E-06,  .869373-05 
. 1 591 2E-04,-. 1 7531 E-05 ,  . 53875E-04 
410  DATA  -. 1 4023E-08,  . 61 762E-07 , -.821 38E-06 ,  .945755-95 
. 1 6593E-04 ,-.1751 0E-05 ,  . 60 1 96E-04 
420  DATA  -.171 35E-08,  .72488E-07 , -.92678E-06 ,  .10247 E-04 
. 1 7266E-04 , - . 1 7489E-05 ,  .61 4835-04 
430  DATA  -.20740E-08,  .84430E-07 , -. 10398E-05 ,  .110603-04 
. 1 7934E-04 , - . 1 7467E-05 ,  . 6273SE-04 
440  DATA  -.24884E-08,  .97653E-07, -. 1 1 604E-05 ,  . 1 1 3955-94 
. 1 8594E-04 , - . 1 7444E-05 ,  • 6 39625-04 
450  DATA  -.2961 6E-08,  . 1 1 222E-06 ,-. 1 2887E-05 ,  .12753E-04 
. 1 9248E-04 , - . 1 7422E-05 ,  . 65 1 535-04 
460  DATA  -.349885-08,  . 1 2820E-06, -. 1 4250E-05 ,  .136323-04 
. 1 9896E-04 ,  - .  1 7399E-05 ,  . 66327E-04 
470  DATA  -.410513-08,  . 1 4566E-06 , -. 1 5691 3-05 ,  .145532-04 
. 20537E-04 , - . 1 7375E-05 ,  . 67470E-04 
430  DATA  -.47859E-08,  . 1 6464E-06 ,-. 1721 25-05 ,  .15452E-04 
.211 71 E-04, -. 1 7353E-05 .  .68538E-04 
490  DATA  -.55464E-08,  . 18522E-06 , -. 1 381 4E-05 ,  .163905-04 
. 2 1 799E-04 ,-.17331 E-05 ,  • 696845-04 
500  DATA  -.63918E-08,  . 20744E-06 , -. 20495E-05 ,  .173475-04 
. 22420E-04 ,-.1731 0E-05 ,  .70757E-04 
510  DATA  -.73276E-08,  . 231 35E-06 , -. 22256E-05 ,  .133215-04 
. 23035E-04 , - . 1 7287S-05 ,  .71 809E-04 
520  DATA  -.83591 E-08,  . 257005-06, -.240985-05,  . 1931 'E-04 
. 23645E-04 , - . 1 7264E-05 ,  . 72O40E-04 
530  DATA  -.9491 5E-08,  . 28444 E-06 ,-.2601 SE-05 ,  .205!  >2-04 
. 24243E-04 ,-.17241 E-05 ,  . 73852E-04 
540  DATA  -.107305-07,  . 31 370E-06,-. 2301 35-05,  .315405-04 
.24344E-04.-.17217E-05,  -74845E-04 
550  DATA  -.130805-07,  .344825-06, -.300965-05,  . 72 4 r. E-04 
• 25435E-04, -. 1 71 92E-05 ,  .7581 9E-04 
560  DATA  -.1 35475-07,  .377865-06, -.333545-05,  .2543V2-04 
. 260205-04 ,-.171 67E-05 ,  • 707775-04 
570  DATA  -.151365-07,  . 41 285E-06,-. 344395-05,  .244375-04 
. 266005-04 ,-.17141 E-05 ,  . 77 ( 1 '  ‘E-04 
580  DATA  -.168525-07,  .44981 3-06 ,-. 36801 5-05 ,  .288635-04 
.271 73E-04,-. 171 155-05,  .786425-04 
590  DATA  -.187005-07,  .438765-06 ,-.39 1395-05,  .286488-04 
. 2774 1 E-04 , - . 1 70885-05 ,  • 79551 5-04 
600  DATA  -.206858-07,  .529735-06, -.41 6525-05,  .277455-04 
. 283025-04 , - . 1 706 1 E-05 ,  . 8 J445E-04 
610  DATA  -.228125-07,  . 57277E-06 ,-. 44 1 908-05 ,  .288535-04 
. 28859E-04 , - . 1 7034E-05 ,  .31 3245-04 
620  DATA  -.250855-07,  .61 7885-06, -.463025-05, 
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.  29409E-04 ,  -  •  1 7006E-05  ,  .  8.?  1 39E-04 
630  DATA  -.27509E-07,  .66508E-06.-.49486E-05,  .310962-04, 
. 29955E-04 , - . 1 6977E-09 ,  .B3040E-04 
640  DATA  - . 30087E-07 ,  . 71 436E-06 , 5224 1 S-05 ,  .322302-34, 
. 30494E-04 , - . 1 6943E-05 ,  .83878E-04 
650  DATA  -.32Q24E-07,  .76576E-06 , -.55066E-05  ,  .  333'/’2K-04  , 
. 3 1 028E-04 ,-.1691 8E-05 ,  • B4702E-04 
660  DATA  -.35727E-07,  .81 932E-06 , --57962E-05 ,  . 34 5 22:: -04 , 
.31 556E-04 , - . 1 6889E-05 ,  .8551 4E-04 
670  DATA  -.3S798E-07,  .87504E-06 , -.60928E-05 ,  .366802-04, 
. 32079E-04 , - . 1 6859E-05 ,  .863 1 4E-04 
680  DATA  -. 42040E-07 ,  . 93289E-06 , -. 63960E-05 ,  .36843.3-04 , 
. 32596E-04 , - . 1 6829E-05 ,  .871 02E-04 
690  DATA  -.45454E-07,  . 99285E-06 , -.67056E-05 ,  .3301 22-04 , 
. 331 08E-04 , - . 1 6799E-05 ,  .87878E-04 
700  DATA  -.49044E-07,  . 10549E-05 , 7021 4E-05 ,  .391352-04, 
. 336 1 4E-04 , - . 1 6768E-05 ,  .08642E-04 
710  DATA  -.52814E-07,  . 1 1 1 91 E-05 , -.73434E-05 ,  .403632-04, 
.3411 6E-04,-. 1 6736E-05 ,  .89395E-04 
720  DATA  -.56767E-07,  . 1 1 855E-05 , -.7671 5E-05 ,  .41544E-04, 
. 346 1 3E-04 , - . 1 6704E-05 ,  -901 37E-04 
730  DATA  -.60903E-07,  . 12539E-05, -.800532-05,  .42729E-04, 
.351 05E-04 , - . 1 6672E-05 ,  . 90868E-04 
740  DATA  -.65224E-07,  .  1 3243E-05  , -.82446E-05 ,  .4391 6E-04 , 
. 35592E-04 , - . 1 6639E-05 ,  .91 589E-04 
750  DATA  -.69730E-07,  .13968E-05,-. 86092 E-05,  .4''106K-04 , 
. 36074E-04 , - . 1 6604E-05 ,  . 92301 E-04 
760  DATA  -.74421 S-07,  . 1 471 2E-05 , -• 90389E-05 ,  .462946-04 , 
. 36550E-04 , - . 1 6568E-05 ,  • 93006E-04 
770  DATA  -.79296E-07,  . 1 5476E-05 , --93934E-05 ,  .474842-34, 
. 37020E-04 , - . 1 6530E-05 ,  . 93705E-04 
780  DATA  -.84357E-07,  . 1 6258E-05 , -.97528E-05 ,  .486742-04, 
. 37483E-04 , - . 1 6490E-05 ,  . 94  397E-04 
790  DATA  -.89603E-07,  . 17059E-05 ,-. 101 1 7E-04 ,  .49362K-04, 
. 3794 1 E-04 , - . 1 6449E-05 ,  • 95082E-04 
800  DATA  -.95031 E-07,  . 1 7878E-05 , -• 1 0485E-04 ,  .51049E-04, 
• 38393E-04 , - . 1 6406E-05 ,  .95761 E-04 
810  DATA  -. 10065E-06,  . 1871 5E-05 , -. 10057E-04 ,  . 52234 E-04 , 
. 33840E-04 ,-.16361 E-05 ,  . 96434E-04 
820  DATA  -.10644E-06,  .  1 9569E-05 ,-•  1 1  233E-04 ,  .57>41  7E-04, 
. 39282E-04 ,-.1631 4E-05 ,  .97101 E-04 
830  DATA  -.11 242E-06 ,  . 20439E-05 , -• 1 1 61 2E-04 ,  .545932-04, 
•  397 1 9E-04 , - . 1 6266E-05 ,  . 97 76 1 E-01 
840  DATA  -.1 1858E-06,  . 21 326E-05 , 1 1 995E-04 ,  .557742-04, 
. 401 52E-04 ,-.1621 6E-05 ,  • 934 1 6E-04 
850  DATA  -.12491E-06,  .22229E-05,-. 12381 E-04,  .569472-04, 
.40578E-04 ,-.161 64E-05 ,  . 900642-04 
860  DATA  -.131 44E-06,  .231 49E-05 1 2771  E-04  ,  .5*1192-34, 
. 4 1 000E-04 ,-.1611 0E-05 ,  . 99707E-04 
870  DATA  -.13814E-06,  . 24034E-05 , -. 1 31 64E-04 ,  .592362-04, 
.4141 8E-04 , - . 1 605 5 E-05 ,  .1 0034 E-03 
880  DATA  -.145013-06,  .25033E-05 , -• 13559E-04 ,  .604492-04, 


. 4 1  a 32E-04 , - • 1 599SE-05 ,  . 1 00975-03 
890  DATA  1 5205E-06 ,  .25996E-05 13956E-04 ,  .6!o05i'-04, 

. 42242E-04 , 1 5939E-05 ,  .101 60E-03 
900  DATA  1 9926E-06 ,  . 26971 E-05 ,-. 14356E-04,  .6275615-04, 

.42649E-04, 1 5879E-08 ,  . 1 0222E-03 
910  DATA  -.16663E-06,  . 27961 E-05 1 4757E-04 ,  .6390215-04, 
. 43032E-04 ,-.1581 7E-05 ,  • 1 0283E-03 
920  DATA  -. 1741 8E-06 ,  . 28964E-05 , -. 1 51 61 E-04 ,  .6504415-04, 
. 43450E-04 , - . 1 5753E-05 ,  . 1 0344E-03 
930  DATA  -.18189E-06,  . 29980E-05  ,  -.  1 5567E-04 ,  .6618015-04, 
•43845E-04, -. 1 568SE-05 ,  . 1 0405E-03 
940  DATA  -. 1 8975E-06 ,  .31 006E-05 , -. 15973E-04,  .6730915-34, 
.44237E-04 , -. 1 5621 E-05 ,  . 1 0464E-03 
950  DATA  -.  1 9776E-06 ,  .32043E-05,-. 16381 E-04,  .6843015-04, 
. 44626E-04 , - . 1 5552E-05 ,  . 1 0524E-03 
960  DATA  -. 20593E-06 ,  .33091 3-05 1 6790E-04 ,  .6954515-04, 
.4501 2E-04 , - . 1 5482E-05 ,  • 1 0582E-03 
970  DATA  -.21423E-06,  . 341 47E-05 , -. 1 7200E-04 ,  .7065215-04, 
. 45396E-04 ,-.1541 0E-05 ,  . 1 0640E-03 
980  DATA  -.22269E-06,  .3521 4E-05  1761 OE-04  ,  .7175515-04, 

. 45777E-04 , - . 1 5337E-05 ,  . 1 0693E-03 
990  DATA  -.23128E-06,  . 36290E-05 , 1 3021 E-04 ,  .72346S-04, 
.461 55E-04 , -. 1 5261 E-05 ,  . 1 0755E-03 
1000  DATA  -. 24001 E-06,  .373738-05,-.  1343215-04,  .739  5015-04, 
.46531 e-04 ,-.151 34E-05 ,  . 1 031 2E-03 
1010  DATA  -.24336E-06,  .38464E-05 , -. 1 3843E-04 ,  .75006E-04, 
. 46906E-04 ,-.151 06E-05 ,  • 1 036SE-03 
1020  DATA  -.2578315-06,  . 39561  E-05 1 9254E-04 ,  .7607215-04, 
. 47279E-04 , -. 1 5027E-05 ,  . 1 0924E-03 
1030  DATA  -. 26692E-06 ,  .4066415-05 1 9663E-04 ,  .7712915-04, 
.47651 E-04, -• 1 4946E-05 ,  . 1 0979E-03 
1040  DATA  -.2761  IE-06,  .41 77215-05, -.20072E-04,  .73W7E-04, 
. 4802 1 E-04 , - • 1 4863E-05 ,  .11 033E-03 
1050  FOR  I  =  1  TO  85 

1060  READ  Cl(l),C2(l),C3(l),C4(l),05(l),C6(l),C7(r) 

1070  NEXT  I 

1080  YK  =  1 

1090  YM  =  YK/1000 

1100  IF  YM  =  10  THEN  YM  =  15 

1110  ’COMPUTE  FA  CLOUT  TIME  OF  ARRIVAL, TA 

1120  TA  =  .1  ’ HOURS 

1130  IF  TA  <  .15  THEN  COSIJB  1510 

1140  ’COMPUTE  g(ta)  [ 1 /Hr] 

1150  GT  =  0 
1160  GOSUB  1270 

1170  L  PR  I  NT  USING  "////.////////  ";TA,GV 
1180  IF  TA  >=  24  THEN  LPRINT  CHR$(l,’) 

1190  IF  TA>-  24  THEN  1230 

1200  IF  TA  <  .1  THEN  TA  =  TA  +  .OlrCOTO  1140 

1210  IF  TA  <  1  !  THEN  TA  =  TA  +  . 1 : GOT > >  1140 

1220  IF  TA  <  24  THEN  TA  =  TA  *  l!:UOT<>  1140 

1230  YK  =  YK  *  10 


A -6 


1240  IF  YK  >  15000  THEN  END  ELSE  1090 
1250  ’ 

1260  ***************** 

1270  ’SUBROUTINE  G(ta) 

1 280  * **************** 

1290  ' 

1300  'COMPUTE  HEIGHT  OF  CLOUD  CENTER, HC  IN  KILOFEET 

1310  HC  »  44 !+6.1*LOG(YM)-.205*(LOG(YM)+2.42)*ABS(LOO(YM)+2.42) 

1320  ’COMPUTE  R  IN  METERS 

1330  HM  =  HC*1 609/5.28  'METERS 

1340  N  =  HM  \  200 

1350  R3  =  Cl (n)/TAa5+C2(N)/TA*4+C3(n)/TA‘3+C4(n)/TA*2+C5(n)/TA 
+C6(N)+C7(N)/SQR(TA) 

1360  RG  =  Cl (N+1 )/TAa5+C2(N+1 )/TA‘4+C3(N+1 )/TAa3+C4(N+1 )/TAa2 
♦C5 (N+1 )/TA+C6(N+1 )+C7(N+1 )/SQR(TA) 

1370  'INTERPOLATE  TO  FIND  R 
1380  R  -  RG-((N+1 )-HM/200)*(RG-RS) 

1390  'COMPUTE  A(r),  AR 
1400  RM  »  R*1E+06 
1410  A2  =  A0+2*B*B 
1420  A3  -  A0+3*3*B 

1430  AR  *  (1 /(SQR(2*PI)*B*RM))*(FV*EXP(-.5*((L0G(RM)-A3)/B)a2) 

♦ ( 1 -FV )*EXP(-. 5*( (LOG (RM)-A2)/B) a2 ) )  ' [ 1 /pm] 

1440  'COMPUTE  dR/dT 

1450  DS  -  -5*C1  (N)/TAa6-4*C2(N)/TAa5-3*C3(N)/TAa4-2*C4(.N)/TAa3 
-C5(N)/TAa2-.5*C7(n)/TAa1 -5 

1460  DG  -  -5*C1 (N+1  )/TAa6-4*C2(N+1 )/TAa5-3*C3(N+1 )/TAa4 

-  2*C4 (N  +  1  )/TAa3-C5(N  +  1 )/TAa2-.5*C7(N  +  1 )/TAa1 .5 
1470  DR  -  d;-((N+1 )-HM/200)*(DG-DS)  ’ [meters/Hr] 

1490  DR  =  DR*1 E+06  '[ pm/Hr ] 

1490  GT  *  AR*ABS(DR)  ' [ 1 /Hr] 

1500  RETURN 
1510  ’ 

1520  ************************ 

1530  'SUBROUTINE  PRINT  HEADER 

1540  ************************ 

1550  ' 

1560  LPRINT"  DATA  FOR  N$; "  PARTICLE  SIZE  DISTRIBUTION" 

1570  LPRINT"  YIELD  =”;YM  "MT" 

1580  LPRINT  ”” 

1590  LPRINT  "  TIME  OF  SMEAR 

1600  LPRINT  "  ARRIVAL  G(T) 

1610  LPRINT  "" 

1620  RETURN 


The  material  in  this  appendix  is  organized  in  four  major  sections, 
one  for  each  particle  size  distribution.  Each  section  contains  a  table 
of  Laurent  coefficients,  two  or  three  tables  of  polynomial  coefficients 
and  five  figures. 

The  Laurent  coefficients  are  used  in  the  function 


g(t)  =  I  ct  t1"1 

i=1 

where  g( t)  is  the  fractional  rate  of  radioactive  fallout 
deposition  per  hour 
t  is  the  time  of  arrival  in  hours 


The  polynomial  coefficients  are  used  in  the  function 


Cj.  -  I  Kj  [  ln(YM)  ]  j-1 


where  C •  is  the  g(t)  Laurent  coefficient 
YM  is  the  nuclear  yield  in  megatons 
n  is  the  degree  of  polynomial 


Each  figure  illustrates  the  g( t)  function  as  computed  by  the  4FIT 
smearing  model,  by  using  eqn  (2)  and  by  using  the  coefficients  for  each 
degree  polynomial  in  eqn  (3). 


r  -  av-  ^ ^ ^ 


Table  B-3 

3rd  Degree  Polynomial  Coefficients 
for  DELFIC  Default  Distribution 


c 

K1 

K2 

K3 

K4 

1 

-.00120380 

-.00001598 

.00000807 

.00000241 

2 

. 1 4890350 

.001 18664 

-.00014724 

.00003443 

3 

.07865119 

-.00207473 

-.00058010 

.00001450 

4 

-.11945574 

-.00953216 

.00151900 

.00016882 

5 

.05029468 

.00642302 

-.00058696 

-.00010676 

6 

-.00935487 

-.00147684 

.00007794 

.00002143 

7 

.00065158 

.00011589 

-.00000296 

-.00000144 

Table  B-4 

2nd  Degree  Polynomial  Coefficients  for 
DELFIC  Default  Distribution 


C 

K1 

K2 

K3 

1 

-.00112340 

.00000003 

- . 00000709 

2 

.15005157 

.00141522 

-.00036363 

3 

.07913485 

-.00197843 

-.OOOO7126 

4 

-.11382593 

-.00841125 

.00045792 

5 

.04673433 

.00571415 

.00008408 

6 

-.00864028 

-.00133456 

-.00005674 

7 

.00060359 

.00010633 
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DELFIC  g(t),  15MT 


*\  *v 

-P  -P 

J  -P  -  P  -P 

oL 

0  -U  T)  “0 
e  -p  0  0  0 
tn  c  j+)*> 

0  P  P  P 
E— 1  L  Q_  CL  CL 
— '  P  e  E  £ 

L_  O  O  O  O 

0  JQOQ 

i  ii  ii  i  n 

□  o  <  +  x 


•0  01  ’0 
( jnoi)  J9d) 


2.0  4.0  6.0  8.0  10. 0  12.0  14.0  16.0  18.0  20.0  22.0 

Time  (hrs) 

Fig.  B-2  DELFIC  g(t),  1MT 
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Fig.  B-5  DELFIC  g(t) ,  1KT 
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Table  B-7 

3rd  Degree  Polynomial  Coefficients 
for  High  Yield  Standard  Distribution 


HI 

K2 

K3 

K4 

0.3  to 

2.0  Hours 

-.12683299 

.01603440 

.00693177 

.  000  ^5 1 0 

1 .05407010 

-.03709439 

-.03059513 

-.001  5'  '905 

-1.29039134 

-.04447673 

.07151674 

.00681682 

•70930171 

.04021040 

-.04955120 

- .00479401 

-.20022912 

-.00351467 

.01663904 

.00146439 

.02836140 

-.00285579 

-.00282655 

-.00017430 

-.0015681 1 

.00053756 

.00019731 

.0000, 144 3 

2.0  to 

24  Hours 

-.00016787 

.00023201 

-.00012187 

- .  000<  >2110 

-.02337136 

-.01052976 

.00639659 

.00003662 

4.70628517 

.63519007 

-.16249243 

-.01998038 

-28.52362613 

-6.47146946 

1  .34574654 

.20414209 

97.72966520 

25.04537544 

-5.39643250 

- . 9052/624 

■176.19628656 

-57 . 28467673 

10.51620436 

1 .87136495 

7  124.88825047  43.75184793  -7.78892336  -1.4412296 6 
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Fig.  B-8  High  Yield  Standard  g(t)  ,  100KT 


Standard  g(t) ,  10KT 


Table  B-10 

3rd  Degree  Polynomial  Coefficients 
for  Low  Yield  Standard  Distribution 


K1  K2  K3 


,00047979 

.00002173 

-.00002986 

,  1 4700891 

.00032719 

.00026261 

,02673082 

-.00390927 

-.00143655 

,08387260 

-.00594650 

.00271652 

>03983654 

.00489934 

-.00134429 

,00788731 

-.00121116 

.00027021 

K4 

.00000208 
.00006540 
.0001 3427 
.00032086 
.0001 9504 
.00004392 


00057176 


00009892  -.00001956  -.00000341 
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Ficj.  B-ll  Low  Yield  Standard  g(t)  ,  15MT 
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Tanle  P-13 

3rd  Degree  Polynomial  Coefficients 
for  NRDL-C61  Distribution 


c 

K1 

K  2 

K3 

K4 

1 

-.00044029 

.00004007 

-.00000866 

-.00000130 

2 

.07706955 

-.00025634 

.00023248 

.0000’ 91  2 

3 

.19412139 

.01671333 

-.00235472 

-.00029320 

4 

-.11804871 

-.02024122 

.00168730 

.00033  31  2 

5 

.03369236 

.00825943 

-.00032405 

-.0001 5130 

6 

-.00472635 

-.00148780 

.00000235 

N-'S 

CO 

6 

o 

q 

r' 

7 

.00025773 

.00009914 

.00000287 

-.000001 38 

Table  P-14 

2nd  Degree  Polynomial  Coeffioi>: 
for  NRDL-C61  Distribution 

■nts 

c 

K1 

K2 

K3 

1 

-.00049035 

.0000301 1 

.000 00078 

2 

.07970812 

.00026900 

o 

1 

3 

.18427774 

.01475334 

-.0004)935 

4 

-.10477230 

-.01759786 

-.0002 1497 

5 

.02864695 

.00723488 

.. 900.,  ..'489 

6 

-.00393167 

-.00132959 

-.0001  47 4  5 

7 

.00021 180 

.00008999 

.<>0001  153 
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Fig.  B-19  NRDL-C61  g(t),  loKT 
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